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Abstract 

We develope a difference calculus analogous to the differential geometry by trans- 
lating the forms and exterior derivatives to similar expressions with difference op- 
erators, and apply the results to fields theory on the lattice [Ref. 1]. Our approach 
has the advantage with respect to other attempts [Ref. 2-6] that the Lorentz in- 
variance is automatically preserved as it can be seen explicitely in the Maxwell, 
Klein-Gordon and Dirac equations on the lattice. 

1 A difference calculus of several independent 
variables 

Given a function of one independent variable the forward and backward differences are 
defined as 

A fix) = fix + Ax) - fix) , Vfix) = fix) - fix - Ax) 
Similarly, we can define the forward and backward promediate operator 

Afix)^^{fix + Ax) + fix)} , Vfix)^^{fix-Ax) + fix)} 

Hence the difference or promediate of the product of two functions follows: 

A{fix)gix)} = Afix)Agix) + Afix)Agix) (1.1) 
A{fix)gix)} = Afix)Agix) + ^Afix)Agix) (1.2) 

This calculus can be enlarged to functions of several independent variables. We use 
the following definitions: 

^xfix.y) = fix + Ax,y) - fix.y) 
^yfix,y) = fix,y + Ay) - fix,y) 

^xfix,y) = ^{fix + Ax,y) + fix,y)} 

\f{x,y) = ]^{fix,y + Ay) + fix.y)} 
^fix,y) = fix + Ax,y + Ay)- fix.y) 
Afix,y) = ^{fix + Ax,y + Ay) + fix,y)} 



These definitons can be easily generalized to more independent variables but for the 
sake of brevity we restrict ourselves to two independent variables. Prom the last definitions 
it can be proved the following identities: 

Af{x,y) = A,Ayf{x,y) + A,Ayf{x,y) (1.3) 

Af{x,y) = A,Ayf{x,y) + ^A,Ayf{x,y) (1.4) 

We can also construct the difference calculus for composite functions. For the sake of 
simplicity we restrict ourselves to functions of two dependent variables and two indepen- 
dent ones, f (u{x,y),v{x,y)). 

We define: 



A„/ = 


f{u + Au,v) - f{u,v) 




A„/ = 


f{u,v + Av) - f{u,v) 




A„/ = 


^{f{u + Au,v) + f{u,v)} 




A„/ = 


^{f(u,v + Av) + f{u,v)} 




A./ = 


f {u{x + Ax, y),v{x + Ax, y)) - f {u{x, y),v{ 




\f ^ 


f {u{x, y + Ay),v{x, y + Ay)) - f {u{x, y),v{ 


x,y)) 


A./ = 


-{/ (u{x + Ax, y),v(x + Ax, y)) + f (u(x, y) 


,v(x,y))} 




^{f {u{x, y + Ay),v{x, y + Ay)) + / {u{x, y), 


,v{x,y))} 



from which the following identities can be proved: 

A/ = AuAJ + AuAJ = A^Ayf + A^Ayf 
Af = AuAJ+^AuAJ = A^Ayf+^A^Ayf 

We can define also the difference operators 

A«,/ = .f(u + A^u,v)-f{u,v) 
^uyf = fiu + AyU,v) - f{u,v) 
^vJ = .fiu,v + A^v)-f\u,v) 

= f{U,V + AyV) - f{U,V) 

and similarly for the promediate operator 

A. J = ^{f{u + A,u,v) + f{u,v)} 
Kyf = ]^{f{u + AyU,v) + f{u,v)] 

= ^{f{u,v + A,v) + f{u,v)} 
Kf = \{f{u,v + Ayv)+f{u,v)} 



From which we deduce the following identities: 



A J = A„,A,J + A„,A„J (1.5) 

AJ = A„^A,J + A„^A,J (1.6) 

AJ = K^A^J + ^A^^A^J (1.7) 

Ayf = A„^ A,J + 1a„^ A, J (1.8) 



These formulas can easily be applied to vector-valued functions: 

M = {ui{x),U2{x), . . .,Un{x)) = u{x) 

and its "tanget vector" 

Au ( Au\ Au-2 Aun 
Ax 



I Ax ' Ax ' ■ ■ ■ ' Ax j ~ 



An inmediate application is the four-position and four-velocity vectors in special rel- 
ativity: 

x'^(r) ^ {x\t),x\t),x'(t),x\t)) 

= /Ax^ Ax^ Ax^ Ax^\ 
~ V Ar ' Ar ' At ' At y 

These vector-valued vector can be expressed as 
for a given set of orthonormal vectors Ca 

2 Discrete differential forms 

Given a vectorial space over Z we can define a real- valued linear function over Z 

f{u) = {uj,u) ueV (2.1) 

The forms u constitue a vectorial hnear space (dual space) *V'^, and can be expanded 
in terms of a basis a;" 

CO — Oa.io'^ 

The basis of and a;° of *V'^ can be contracted in the following way 

(a;", 6^3) 5^ (2.2) 

hence 

(a;,ea) = (Ta , ,u) = , {uj,u) = aau'^ (2.3) 



If we take ujf^ = Ax'^ as coordinate basis for the linear forms we can construct discrete 
differential forms (a discrete version of the continuous differential forms) 

A particular example of this discrete form is the total difference operator (1,3) of a 
function of several discrete variables written in the following way: 

A/(^, y) - (^) A. + Ay (2.4) 

For these discrete forms we can define the exterior product of two form a and S 

p Aa — —a A p 

which is linear in both arguments. 

For the coordinate basis we also have 

Ax A Ay = -Ay A Ax 

With the help of this exterior product we can construct a second order discrete differ- 
ential form or 2-form, namely 

p A (7 = -paAx" A Uo^Ax^ = \- {pocOp - ppUa) Ax"' A Ax^ = GocpAx"' A Ax^ (2.5) 

where Oap is an antisymmetric tensor. Similarly we can define a discrete p-form in a 
n-dimensional space (p < n) 

^ = Yfni2 ■ ■ ■ iA^'' A Ax'^ ... A Ax'- 

where Ci^i^ ■ ■ - ip is a totally antysymmctric tensor 

The dual of a p-form in a n-dimensional space is the (n — p)-form *a with components 

— iLn''i*2-«P<r . . . , , 

where e is the n-dimensional Levy-Civitta totally antisymmetric tensor (£123... = 1) 
We give now some examples: 

Energy-momentum 1-form 

P = -EAt + P^Ax + PyAy + P^Az (2.6) 

where (P^, Py, Pz, iE) = Pn is the four-momentum. 

Vector potential 1-form 

A = Af.Ax'^ = Aa^Ax + Ay Ay + A^Az + At At (2.7) 

where — {A^, Ay, Az, At) is the four-potential. 



Charge- current 1-form 

J = J^Ax'' = Jj^Ax + Jyi^y + J^/^z - pAt (2.8) 
where is the density current four-vector. 
Faraday 2-form 

F = E^Ax A At + EyAy A At + E^Az A At 

1 

+ B^Ay AAz + ByAz A Ax + B^Ax A Ay = -F^^Ax'' A Ax" (2.9) 

with {Bx, By, B^) = B and {E^, Ey, E^) = E the magnetic and electric field, respectively. 
Maxwell 2-form (dual of Faraday 2-form) 

*F = ^^s^^xkF^'^^x'' a Ax"" = - B^Ax A At - ByAy A At - B^Az A At 

E^AyAAz^EyAzAAx^E^AxAAy{2.lQi) 

3 Exterior calculus 

Given a 1-form in a two-dimensional space 

uj — a{x, y)Ax -\- b{x, y)Ay 
we can define the exterior difference, in the similar way as the exterior derivative, namely, 
Aa; = AaAAx + AhA Ay 

\ Ax Ay j \ Ax Ay ) 

= \^ir^ \ Ax A Ay 3.1 

V Ax Ay y ^ ^ > 

where in the last expression we have used the properties of the exterior product. 

This definition of exterior difference can be easily written for 1-form in n-dimensional 

space. 

Given a 2-form in a 3-dimensional space, 

u = a{x, y, z)Ay A Az-\- b{x, y, z)Az A Ax + c{x, y, z)Ax A Ay (3.2) 
we can also define the exterior difference as: 

Aa; = AaAAyAAz + AbAAzAAx + AcAAx A Ay 

( A^AyA^a A,ASJ) A^A,,A^c\ , , . , , 

= / + " " + / Ax A Ay A Az 3.3 

I Ax Ay Az I ^ ^ ^ 



Given a 3-form in a 4- dimensional space 

uj = aAy A + + hAz A Ai A Ax + cAt AAxM^y + dAx A Ay A Az 

we can define an exterior difference as before: 

cu = AaA Ay A Az A At + A6 A A2; A At A Ax 

+Ac A At A Ax A Ay + Ad A Ax A Ay A Az 
f A^AyA.Ata A^AyA.Atb ^ A^AyA.Atc 



Ax Ay Az 

Ax A Ay A Az A At (3.4) 



A^A^A.A^d 



At 

Tfie exterior derivative applied to the product of a 0-form (scalar function f) and a 
1-form (a; = aAx + 6Ay) is 

A {fu) = AfAu + Af AAu (3.5) 

where A/ is expressed in (1.4) and Auo = AaAx + A&Ay 

The exterior difference of the product of two 1-forms is easily obtained 

A {coi A UJ2} = Aci;i A Auj2 - Aui A Alo^ (3.6) 
The exterior difference of the product of a p-form p and a q-form a is 

A{pA(7} = ApAA(7 + (-l)^ApAA(7 (3.7) 
Finally for any p-form uj we have 

A^uj = A (Aa;) = (3.8) 

Some examples: 

Prom the Faraday 2-form we write down one set of Maxwell difference equations 

AF = A (AA) = 



f A^AyA^AtB^ A^AyA^AtBy A^AyA^AtB^\ , 

\^ * + \" * ^ + y ^ * ' Ax A Ay A Az 

V Ax Ay Az y 

/ A,A,A,A,E, A,A,A,A,£;, A,A,A,A,ig, \ ^^^.^^^ 
\ At Ay Az J ^ 

fA^A^A^AtBy A^AyA^AtE^ A^AyA^AtEA , 

X y z t y X y z t x x y z t z W ^ A Ax 

V At Az A.T ; 

(A^AyA-^AtB^ Aj^A^A^AjE^ A^A„A^At£^^\ , , , , , 

+ ^ y t ^ ^ X y z t y X y z ^ At A Ax A Ay 3.9 

I At Ax Ay y \ ) 



from the Maxwell 2-dual form *F we get the other set of Maxwell equations: 

A*F = 47r* J 

\ Ax Ay Az j 

(A^A^A AtEx AxAyA^tBz A^AyA^AtBA , 

X y z t X X y z t z x y z t y At A Aw A A^ 

\ At Ay Az J ^ 

( AxAyA^AtEy AxAyA^A.Bx A,A,AA^, \ 

At Az Ax J 

fA^AyA AtEz A^AyA^AtBy A^AyA.AtBA , 

X y z t z X y z t y x y z t x At A Ax A Aw 

\ At Ax Ay j ^ 

= 47r {pAx A Ay A Az - J^At A Ay A Az 

- J^^At A Az A Ax - J;,At A Ax A Ay) (3.10) 

Taking the exterior derivative of the last equation we get an other example of A^ = 0. 

(AxAyA^Atp A^AyA^AtJx A^AyA^AtJy A^^A^A^A^JA ^ ^ ^ ^ 

^ + + ^ + -AtAAxAAyAAz = 

V At Ax Ay Az I ^ 

(3.11) 

Note that the coefficient of the difference form is the discrete version of the continuity 
equation. 

From a scalar function we get the wave equations in terms of difference operators, 
namely, 

-*A*A(\) = U(^ (3.12) 
where □ is the discrete d'Alambertian operator: 

{-V.V.V.Vt (A,A,A,Ai) + V Vt [A.,AyA,At) 

+VxVyV,Vt (Ax^yA.At) + VxVyV.Vt [a^ Ay A, At) } (xyzt) = (3. 13) 

From the vector potential 1-forms A = A^Ax^^ we can construct Faraday 2-form 
F = AA from which the Maxwell equations are derived 

A*F = A*AA = 47r*J 

Taking the dual of this expression we obtain 

*A*AA = 47rJ (3.14) 

If we choose the Lorentz condition 

A^AyA^A,^^ ^ A^A^A,A,Ay ^ A^A^A.At^, _ A^A^A^A^A ^ q .3 

Ax Ay Az At ^ ' ^ 

we obtain the wave equation for the vector potential 

UA^^AttJ^ (3.16) 
where □ is the d'Alambertian defined in (3.13) 



4 Lorentz transformations 

In order to compute the transformation of the discrete differential forms we start with 
the coordinate-independent nature of 1-form 

oj = Uf.Ax" (4.1) 
where the Ax'^ are the space-time intervals in Minskowski space-time. From 

Ax^' = A'^'Ax'' (4.2) 

where A^' is a global Lorentz transformation, and from the coordinate-free expresion for 
cu we get 



co^^f = c^.A;, (4.3) 
From the total difference of a function of several variables / {x, y, z, t) 



Recall that A^'.A'^' = S'! 



^ A^AA/ Ax + h^y^^/^y + M^AA/ + MAM At (4.4) 
Ax Ay Az At 



it follows that the coefficients of the 1-forms, namely, 

A, A A A A,\ 

(4.5) 



'A^AyA.Atf A^AyA.Atf A^A^A.A^ A^A^A^A^^ 



\ Ax Ay Az At J 

transform covariantly like the coefficients a;^ of (4.1). Note that in this case the meaning 
of Ax, Ay, Az, At in the denominator is different from the meaning in the numerator, 
because the later Ax^ are elements of the exterior products, and the former Ax are small 
scalars. The different roll of these quantities becomes clear in the continuous limit, where 

Ax ^ ds and ^2%^ ^ f 

Ax ox 

For the 2-form in Minkowski space 

F = ^F^^Ax^ A Ax" 

we obtain 

F = ^F^^.Ax'" A Ax"' 

because of the coordinate independent nature of the Faraday 2-form. From the transfor- 
mation of Ax'^ (see (4.2)) and the properties of the exterior product we get 

F^,., - F^pA'^^Ai (4.6) 

The same technic can be applied to components of discrete p-forms. The dual of a 
p-form are also discrete (n — p)-form, therefore their components transform covariantly 
like totally antisymmetric tensor. For instance the components of Maxwell 2-forms. 



{*F)^^ = e^p,^F^^ (ei234 = +1) 

transform covariantly 

CF)a,0, = CF).xKA% (4.7) 

With this definition of the transformation of the components of a discrete p-form we 
can prove the covariance of the discrete wave equation, and the covariance of Maxwell 
equation in discrete form 

5 Application to Klein-Gordon and Dirac wave 
equation on the lattice 

We define the scalar function on the (3+1) dimensional cubic lattice 

(jid, i2e2, ises, nr) = (f){j,n) 

where ei,e2,e3,r are small quantities in the space-time directions and ji,j2,j3,'n are in- 
teger numbers. 

We define the difference operators 

Then the Klein- Gordon wave equations defined on the grid points of the lattice can 
be read off 

(5+5^ + 5+6^ + 5+5^ - St 5^ - M^ri+ri~) (j, n) = (5.1) 

It can be verified by direct substitution that the plane wave solution satisfy the dif- 
ference equation 

"^'"J- l,l-ii£,fcj [l-^teoM] [l-ile,k,) [l + iiTu) '■"■■'> 

provided the "dispersion relation" is satisfied 

cu^ -kf-kl-kl^ (5.3) 



From the last section the Klein- Gordon equation is invariant under finite Lorentz 
transformations. Imposing boundary conditions on the plan waves we can construct a 
complete set of orthogonal functions, hence a Fourier analysis can be developped as it has 
been done in Ref. [1]. 

The discrete version of the Dirac wave equation can be written as 

(71^1+ + l2^t + 7353"" - ^lA^t + ^^^) (j, = (5.4) 
where 7^, = 1, 2, 3, 4 are the usual Dirac matrices. Applying the operator 

(liST + 12^2 + 73^i" - ilA^I - Mr]~^ 

from the left on both sides of (5.4) we recover the Klein-Gordon equation (5.1). Let now 
construct solutions to (5.4) of the form 

il; {j, n) = LV (k, E) f {j, n) 

where the / {j, n) are given in (5.2). 

The four-component spinors u (^k, E^, with momentum k = {ki,k2,k3), must satisfy 

(i^-k- -f^E + m) uo (k, e) =0 (5.5) 
as in the continuous case. Multiplying this equation from the left by 

(i^-k- 74^ - m) 
we obtain the dispersion relation 

E^-k^^ (5.6) 

From (5.5) we obtain the spinors solutions corresponding to positive energy as in the 
continuous case 

- (l w\ Et + M ^l^ ( \ K 7^ 



with Ek = + V A;^ + M2, (7=1,2 and 
Similarly for negative energy spinors. 

In order to analyze the Lorentz invariant of the Dirac wave equation on the lattice we 
take the Lorentz transformations of the difference operators as in (4.3) 



= (5.8) 
5;' = A,j;- (5.9) 

The linear transformation law for tp 



is determined by the requirement that ipi satisfy the same equation in the transformed 
frame as does ^l) in the original frame 

(7;.5+' + m)V'/ = (5.10) 

With (5.8) we have 

(A^.7.C + 5V = (5.11) 

Multiplying from the left by S"^ we recover the Dirac equation in the original form 
provided 

5-^^^ = A^,7. (5.12) 

A particular solution for the Lorentz transformations of the spinors in the case of 
rotations is 

■0/ = S''0 = (cos q; + i sin 075) (5.13) 

where a is the angle of rotation and 75 = 71727374. In the case of pure Lorentz transfor- 
mation we have 

4)f = Sip = chPtp + shp-f4'y2'ip (5.14) 
with ip — and (3 is the usual parameter such that thp — ^ 
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